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I. INTRODUCTION 



Liquid crystals (LCs) are states of matter whose properties are in between those of a 
crystaUine sohd and an isotropic hquid phase pQ. They are usually classified in terms of 
positional and orientational order. Nematic LCs exhibit long-range orientational order, 
as the anisotropic particles are on average aligned along a preferred direction, but they 
lack long-range positional order. Smectic phases consist of stacks of fluid-like layers of 
orientationally ordered particles, where each layer is often considered to be a two-dimensional 
fluid. Onsager showed in his seminal contribution the existence of a purely entropy-driven 
isotropic-to-nematic (I-N) transition in a system of infinitely long hard rods [2J. Moreover, 
the I-N transition was confirmed by computer simulations for systems of hard rods with 
finite length. Additionally, Frenkel and coworkers explored the formation of smectic LCs 
of perfectly aligned [3] and freely rotating ^ hard rods, and found a thermodynamically 
stable smectic phase of hard rods as a result of entropic effects. The equilibrium properties 
of smectic LCs are well-studied and are well-understood by now [3]. Experimental P-[H], 
theoretical [HHIl] , and computational fTEj [TB] studies have analyzed the phase behavior and 
structure of smectic LCs of colloidal hard rods. Other investigations involve extensions to 
binary mixtures with rods of different geometry [T7H23] . with other anisotropic [2H - I27] or 
spherical colloidal particles [28143T] . or with non-adsorbing polymer as depletants [32143^ . 

By contrast, the dynamics on a single-particle level in smectic LCs have only recently 
received attention, although an early study on the diffusion (or "permeation") of anisotropic 
particles had already been reported more than forty years ago by Helfrich in order to ex- 
plain the capillary flow in cholesteric and smectic LCs [35j. Substantial advances in new 
experimental techniques (e.g. NMR coupled to strong magnetic field gradients |36j or flu- 
orescent labeling [3YJ) disclosed the non-Gaussian nature and quasi-quantized behavior of 
the layer-to-layer diffusion. These achievements sparked off new theoretical work, based 
on dynamic density functional theory, which not only confirmed the non-Gaussian layer- 
to-layer hopping-type diffusion and the presence of permanent barriers due to the static 
smectic background, but also showed the relevance of temporary cages due to the mutual 
trapping of neighboring particles [381 EH]- We note that non-Gaussian dynamics due to a 
rattling- and-jumping diffusive behavior is common in two-dimensional liquids [40j, cluster 
crystals jH], and glasses [12], and has also been observed for the diffusion of a single par- 
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tide in a periodic external potential [13]. It is therefore not surprising to observe similar 
behavior in smectic LCs. Our simulations of parallel [H] and freely rotating hard rods 
supported indeed these conclusions, but unveiled in addition a striking analogy with the 
non-exponential structural relaxation and non-Gaussian dynamics of supercooled liquids. 
The non-exponential relaxation of the density fluctuations might be due to either a het- 
erogeneous scenario with particles relaxing exponentially at different relaxation rates, or a 
homogeneous scenario with particles relaxing non-exponentially at very similar rates jl6]. 
Here, we investigate the collective motion of fast-moving rods in stringlike clusters. Coop- 
erative diffusion, which accounts for the non-exponential decay of the correlation functions, 
yields an intriguing link between the dynamics observed in equilibrium smectic LCs and that 
of out- of- equilibrium supercooled liquids. As far as smectic LCs are concerned, this remark- 
able collective motion could not be captured by the one-particle analysis of Ref. |13], and 
was not observed in Ref. |17] . Two-dimensional liquids of soft disks [ID] and cluster crystals 
[H] did not show this feature either. By contrast, several experimental and computational 
studies on glassy systems reported the existence of structural heterogeneities [15V[5^ . In 
particular, Glotzer and coworkers performed molecular dynamics simulations on a fragile 
glass-forming liquid and detected cooperative motion of stringlike clusters with an increas- 
ing string length of up to ~ 15 particles by cooling the system towards the glass transition 
|48j . Similar results were observed more recently in silica, a strong glass-former, suggesting 
that stringlike motion is a universal property of supercooled liquids [32] ■ 

In this paper, we investigate the effect of anisotropy of the rods on the non-Gaussian 
layer-to-layer diffusion and cooperative motion of stringlike clusters in bulk smectic LCs of 
freely rotating hard rods. In supercooled liquids, it is generally accepted that in the case of 
heterogeneous dynamics, the cooperative motion of particle clusters plays a crucial role in 
the structural relaxation [53]. We argue that a similar behavior can be observed in smectic 
LCs, where cooperative layer-to-layer motion of strings with various sizes contributes to the 
long-time relaxation behavior of the smectic phase. 

This paper is organized as follows. In section II, we introduce the model, the simulation 
details, and the computational tools to describe the layer-to-layer diffusion. The results on 
the non-Gaussian and heterogeneous dynamics, as well as evidence of cooperative motion 
are discussed in section III. In the last section, we present our conclusions. 
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II. MODEL AND SIMULATION METHODOLOGY 



We perform simulations of a system with = 1530 — 3000 freely rotating hard sphe- 
rocylindrical rods with diameter D and length L + D, distributed over 5-10 smectic layers 
of approximately 300 rods each. Three different aspect ratios, L* = L/D, are considered: 
L* = 3.4,3.8, and 5.0. The region of stability of the smectic phase decreases with L*, and 
disappears at L* < 3.1, where only a stable isotropic-crystal phase transition is found [T4] . 
For L* = 3.8 and 5.0, the smectic phase is stable for 2.3 < P* < 2.8 and 1.4 < P* < 2.3, 
respectively, where P* = PPD^ is the reduced pressure. P is l/Zc^T, with ks the Boltz- 
mann's constant and T the temperature. For lower P*, the smectic phase transforms into a 
nematic phase, while for higher P* the smectic freezes into a crystal phase. For L* = 3.4, 
the nematic phase is unstable, and the smectic phase transforms directly into an isotropic 
phase for P* < 2.8 and crystallizes for P* > 3.0. We studied the dynamics of the bulk 
smectic phase at the pressures and packing fractions indicated in Table I. For convenience, 
we label the systems with Si — Sq (see Table I). 

We performed Monte Carlo (MC) simulations in a rectangular box of volume V with 
periodic boundary conditions. Firstly, we performed equilibration runs in the isobaric- 
isothermal (NPT) ensemble to expand the system from an ordered crystalline phase to an 
equilibrated smectic phase. Each MC cycle consisted of N attempts to displace and/or rotate 
the randomly selected particles, plus an attempt to change the box volume by modifying 
the three box lengths independently. Translational and rotational moves were accepted if 
no overlap was detected. The systems were considered to be equilibrated when the packing 
fraction reached a constant value within the statistical fluctuations. A typical equilibration 
run took roughly 3 x 10^ MC cycles and was followed by a production run in the isochoric- 
isothermal (NVT) ensemble to analyze the relaxation dynamics. At this stage, we kept 
the volume constant to avoid unphysical collective moves which do not mimic the Brownian 
dynamics of the rods properly. Standard MC simulations with small displacements were used 
to mimic Brownian motion. This computational approach was shown to be very efficient to 
study the slow relaxation of glasses at low temperatures [54] or at high concentrations [55] . 
We fixed the maximum displacement according to (i) a reasonable time of simulation, (ii) 
a satisfactory acceptance rate, and [iii) a suitable description of the Brownian motion of 
colloidal particles suspended in a fiuid (see Fig. 1). To this end, we monitored the mean- 
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L* 


3.4 


3.8 


5.0 


Stable Sm 


2.8 < P* < 3.0 


2.3 <P*< 2.8 


1.4 <P*< 2.3 


p* 


2.85 


3.00 


2.35 


2.50 


1.60 


2.00 


V 


0.556 


0.568 


0.536 


0.551 


0.508 


0.557 


h/{L + D) 


1.018 


1.014 


1.023 


1.015 


1.048 


1.030 


ct/{L + D) 


0.050 


0.043 


0.075 


0.052 


0.093 


0.043 




0.036 


0.020 


0.116 


0.050 


0.325 


0.028 




0.230 


0.179 


0.295 


0.256 


0.413 


0.268 


t/T 


0.19 


0.23 


0.17 


0.20 


0.14 


0.20 




0.28 


0.30 


0.27 


0.32 


0.20 


0.27 


J 1 


1.10 


1.35 


1.40 


1.65 


1.00 


1.20 


I 


0.04 


0.02 


0.04 


0.04 


0.09 


0.03 


fc 


0.04 


0.02 


0.05 


0.05 


0.09 


0.03 


jmax 


0.06 


0.04 


0.14 


0.07 


0.25 


0.04 


Uo/kBT 


6.2 


7.5 


4.0 


5.8 


3.5 


7.5 




Si 


^2 


Ss 


5*4 


^5 


5*6 



TABLE I: Details of the systems that we studied in this paper, consisting of hard spherocylinders 
with varying length-to-diameter ratio L* = L/D and reduced pressures P* = (3PD^, and corre- 
sponding packing fractions ry. For comparison, we give the pressure range of the stable smectic 
phase for the corresponding systems. Additionally, we give the layer spacing (h); the standard 
deviation of the displacement from the equilibrium smectic phase (a) in units of (L+D); the long- 
time in-layer D^y and inter-layer diffusion coefficients in units of t/D"^; the most probable (i), 
median (tj), and maximal (t™"^) jump times; the fraction of collective jumps f^, /*, and f^'^^^ 
calculated with a temporal interval Ato equal to t, t*j, and tj*"^, respectively; and finally, the height 
of the energy barriers (Uq) in units of fc^T. The systems are labeled by Si — Sq. 

square displacement in the z and xy directions for several values of the maximum step size 
^max, with 6max,z = '^^max,xy to take into account the anisotropy of the self-diffusion of the 
rods |56j. We found Smax,xy = -D/10 and 6max,z = D/5 to be the optimal values which 
satisfied the above requirements. 
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FIG. 1: Trajectories in a plane perpendicular to the nematic director of approximately 300 rods 
in a smectic layer of hard spherocylinders with a length-to-diameter ratio L* = 5.0 and reduced 
pressure P* = 1.60 collected over 5 x 10^ MC cycles. 

As unit of time, we have chosen r = D'^/D^^., where D^^. is the short-time translational 
diffusion coefficient at L* = 5.0, which is the isotropic average of the diffusion coefficients in 
the three space dimensions: Df^. = {Dl + 2D'^y)/3. At short times, when the single particle is 
rattling around its original position without feeling the presence of its surrounding neighbors, 
the dependence of Df^ on the pressure can be safely neglected. We checked that our results 
(measured in units of r) were independent of S^ax- 

In order to characterize the layer-to-layer hopping-type diffusion and the structural relax- 
ation of our systems, we calculated (z) the energy barrier, {ii) the self part of the van Hove 
correlation function, {in) the non-Gaussian parameter, (iv) the mean square displacement, 
(v) the intermediate scattering function, {vi) the probability distribution of the size of the 
stringlike clusters, and (vii) their dynamic cooperativity. 

Energy barrier. We computed the energy barriers from the (relative) probability 7r(z) of 
finding a rod at a given position z along the nematic director n. As reported in Ref. |57] . 
this probability is proportional to the Boltzmann factor 

n{z) oc exp [~U{z)/kBT] , (1) 
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where U{z) denotes the effective potential for the layer-to-layer diffusion. 

Self-part of the van Hove correlation function. To quantify the heterogeneous dynamics 
due to the rattling and hopping type 2;-diffusion of the rods, we calculate the self-part of the 
van Hove correlation function (VHF), which measures the probability distribution for the 
displacements of the rods along the nematic director h at time to + 1, given their z positions 
at to- It is defined as [57] 



N 



i=l 



with (...) the ensemble average over all particles and initial time to, and S is the Dirac-delta. 
Note that Gs{z,t) would be a Gaussian distribution of z for freely diffusive particles. 

Non- Gaussian parameter. A quantitative description of the non-Gaussian behavior of the 
layer-to-layer diffusion can be obtained in terms of the non-Gaussian parameter (NGP) [58] : 



:i + 2/d) {Az{tyy 

where Az{t) = z{to + t) — z{tQ) denotes the z-displacement of a rod in the time interval t 
starting at to 5 and d = 1. Heterogeneous dynamics occurs on a time scale t, if the NGP is 
non- vanishing. For the in-layer diffusion, a similar NGP, a2,xy{t), with d = 2 can be defined. 

Intermediate scattering function. The structural relaxation and, in particular, the decay 
of the density fluctuations can be quantified by the self-part of the intermediate scattering 
function (ISF) 

F,(t) = (exp[2q-Ar(t)]) (4) 

at the wave vectors q = (0, 0, g^) and {q^, qy, 0) corresponding to the main peaks of the static 
structure factor in the z and xy directions, respectively. 

Cluster size distribution. In order to investigate whether or not there is cooperative 
motion of stringlike clusters of particles, which might be responsible for the heterogeneous 
dynamics and stretched-exponential decay of the ISF, we first determine the clusters of 
fast-moving particles. The fast-moving particles can be identified as those particles that 
have traveled a substantially longer distance than the average in a certain time interval and 
are intimately related to the particles that reside in between the center-of-masses of two 
smectic layers. Hence, particles that reside more than some rattling distance 6rat from the 



center-of-mass of the nearest smectic layer in a static configuration are defined as interlayer 
rods. To define 5rat, we first calculate the variance of one period of 7r(2;), defined as 

rh/2 

= z^7r{z)dz, (5) 

J-h/2 

where h is the layer spacing obtained from fits to the density profiles Tr{z), and '^{^)d^ — 
1. We assume that two interlayer rods belong to the same string if their z and xy distances 
are smaller than h and D, respectively. We calculated the probability distribution P[n) of 
the number of rods n in a string for Srat — 2cr. In Table I, we give the values of the layer 
spacing h and the standard deviation (square root of the variance) a for the systems that 
we studied. We find that denser states reveal a slightly smaller layer spacing for all aspect 
ratios, as expected. 

Dynamic cooperativity. To check whether or not the rods in a given string move col- 
lectively from one layer to another layer, we require a cluster criterion involving a spatial 
proximity of particles and a temporal proximity of jumping rods. To this end, we assume 
that two jumping rods i and j are moving cooperatively if (1) their arrival times t^^^ and 
t^") in their new layers (i.e. the time at which their distance to the middle of the new layer 
equals 6rat) satisfy \t^'^^ — t'^^^ \ < Ato, and (2) the center-of-mass positions of the two particles, 
rj(t'^*^) and rj{t^^^), satisfy the above mentioned static cluster criterion. The time interval 
Ato is fixed on the basis of the distribution of jump times and will be introduced in the 
following section. 

III. RESULTS 

The layered structure of smectic LCs yields an effective periodic potential U (z) for the 
diffusion of rods out of the middle of a smectic layer to another layer. The permanent 
energy barriers for the layer-to-layer diffusion are determined from the effective potential 
U{z) introduced in Eq. (1). In order to evaluate the effect of the pressure on the inter-layer 
diffusion, two separate state points are considered for each aspect ratio: one in the proximity 
of the nematic-smectic (N-Sm) or isotropic-smectic (I-Sm) transition, and the other in the 
proximity of the smectic-crystal (Sm-K) transition. In Fig. 2, we show U{z) for the six 
systems along with the fit U{z) = [sin(7r2;//i)]^*, where m is an integer number. 
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Uo = Yl^i is barrier height, and h the interlayer spacing given in Table I. 




FIG. 2: Potential energy barriers for the layer-to-layer diffusion in smectic LCs of hard spherocylin- 
ders with varying length-to-diameter ratios L* and reduced pressures P* . The open symbols refer 
to weak smectic states L* = 3.4; P* = 2.85 (a), L* = 3.8; P* = 2.35 (□), and L* = 5.0; P* = 1.60 
(O). The solid symbols denote state points deep in the smectic phase L* = 3.4; P* = 3.00 (a), 
L* = 3.8; P* = 2.50 (■), and L* = 5.0; P* = 2.00 (•). The solid lines are fits with m = 5 harmonic 
modes. 

As a general consideration, we observe that the energy barriers increase with increasing 
packing fraction. The denser state is characterized by a much stronger confinement of the 
particles to the middle of the smectic layers. This is especially evident at L* = 5.0 as the 
height of the barrier increases from S.bksT to l.'bksT from 7] = 0.508 to 0.557. A similar 
behavior was detected in experiments on smectic LCs of fd viruses, where the height of 
the energy barriers was found to increase from 0.66A;bT to l.SGksT by decreasing the ionic 
strength ^j. According to the authors, a low ionic strength gives rise to stronger correlations 
due to more pronounced electrostatic interactions between the virus particles, resulting in 
higher energy barriers. Older experiments on thermotropic liquid crystals estimated energy 
barriers of ^l-AksT [SSI ED]- Our results are in good quantitative agreement with those 
obtained by computer simulations in Ref. |6T], where rods with L* = 3.8 and 5.0 were 
studied. Additionally, these authors found that the energy barrier for a rod to achieve 
a transverse inter-layer position was several kBT higher than that to diffuse from layer 
to layer, confirming the difficulty to detect transverse particles in between smectic layers 
[52] . We further observe that the packing fraction is not the only parameter affecting the 
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effective potentials U{z) as displayed in Fig. 2. It is interesting to note that the two 
systems 5*1 and 5*6, with L* = 3.4 and P* = 2.85 as denoted by the open triangles (A) 
and with L* = 5.0 and P* = 2.00 as shown by the solid circles (•), corresponding to 
very similar packing fractions rj = 0.556 and 0.557, respectively, exhibit barrier heights 
that differ by approximately 1.5A;bT. The barrier height at fixed rj thus increases with 
the particle anisotropy. The effect of the particle anisotropy on the effective potential can 
also be illustrated by comparing the barriers for the systems 5*2 and 5*6 with L* = 3.4 and 
P* = 3.00 (a) and L* = 5.0 and P* = 2.00 (•), respectively. In this case, although the 
packing fractions are significantly different [r] = 0.568 and 0.557), both barrier heights are 
7.5kBT. This result shows that a more pronounced particle anisotropy yields significantly 
higher energy barriers. In conclusion, the barrier height increases with increasing packing 
fraction and particle anisotropy. Moreover, in our recent work on perfectly aligned hard rods, 
we noticed that the freezing out of the rotational degrees of freedom has also a quite tangible 
impact on the height of the barriers, which were found to be higher than those observed in 
systems of freely rotating rods, especially at low packing fractions [H]. Such barriers fade 
out gradually by approaching the continuous N-Sm transition, while they remain finite in 
the case of the first order N-Sm transition of freely rotating hard rods. One might expect 
that structural defects, such as screw dislocations and stacking faults, could facilitate the 
layer-to-layer diffusion by creating barrier-free nematic-like pathways through the layers |63] . 
The small size of our system and the periodic boundary conditions preclude the development 
of such defects. 

The periodic shape of the effective potential determines a hopping-type diffusion along 
the direction parallel to the nematic director n, with the particles rattling around in their 
original layer until they find the appropriate conditions to overcome the barrier and jump 
to a neighboring smectic layer. An efficacious way to quantify the diffusion of rattling and 
jumping rods along n is provided by the computation of the self part of the VHFs defined in 
Eq. (2). In Fig. 3, we show the self VHFs for the six systems of interest as a function of z 
at several equidistant times t. We detect the appearance of peaks at well-defined locations 
corresponding to the center-of-mass of the smectic layers along n, in agreement with previous 
experimental [37] and theoretical [38] results. For each aspect ratio, we note that the height 
of the peaks is larger and the number of peaks is generally higher at the lowest packing 
fraction. This indicates that decreasing the pressure leads to a higher number of jumping 
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particles which are able to diffuse longer distances. These fast particles determine the 
heterogeneous dynamics of the system and affect its structural relaxation. 

10-^ 
10-^ 

^ 10"^ 

10-^ 
10-^ 

-2-1012 -2-1012 

z/h 

FIG. 3: Self part of the van Hove function Gs{z,t) for smectic LCs of hard spherocylinders with 
varying length-to-diameter ratios L* and reduced pressures P* . Length-to-diameter ratio, pressure, 
and packing fraction are given in each frame as {L* , P* , rf). The curves refer to the time evolution, 
from t = 0.4r (dotted lines) to t = 40r (solid lines) with increments of w Sr. 

Deviations from the Gaussian behavior of the VHFs have been extensively analyzed in 
liquid [lOl ini] , glassy [l2l HS] , and liquid crystalline [381 El] systems in terms of the non- 
Gaussian parameter defined in Eq. (3). The NGPs as measured in our six systems are 
shown in Fig. 4. The in-layer NGPs, a2^xy, are essentially negligible for the whole time 
range, implying a Gaussian in-layer dynamics which is typical of a liquid-like system. The 
layer-to-layer NGPs exhibit a time-dependent behavior, which is strictly linked to the caging 
effect exerted on the rods by their nearest neighbors. More specifically, 02,2 is basically 
zero at short times when the rods are rattling around their original location and do not 
perceive the presence of the surrounding cage. Between t/r = 0.1 and 1.0, ^2,2 starts to 
increase, indicating the development of dynamical heterogeneities. During this time interval, 
the motion of the rods is hampered by the trapping cages and becomes subdiffusive. It is 
reasonable to assume that the average life-time of the cages corresponds to the time between 
0:2,2 — and 0^2,2 = Oi^^T 1 t^e maximum value of the NGP, which occurs at time t = tmax with 
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2 < tjnax/T ^ 10, depending on the system. This peak increases with packing fraction and 
its location determines the beginning of the long-time diffusive regime, where the deviations 
from Gaussian behavior start to decrease. Increasing the pressure affects the z-diffusion in 
a twofold manner: («) it increases its heterogeneous behavior, and {ii) it delays the onset of 
the long-time diffusive behavior. Comparison of the two systems 5*1 and 5*6, corresponding 
to very similar packing fractions, shows that the dynamics becomes more non-Gaussian for 
system Sq consisting of longer rods and higher energy barriers. On the other hand, for 
systems S2 and Sq, displaying similar barrier heights, the non-Gaussian dynamics is again 
more pronounced for system 5*6 with longer particles. Hence, increasing the anisotropy of 
the particles yields higher energy barriers and dynamics that is more heterogeneous and 
non- Gaussian. 

6 

3 


6 

S 3 


6 

3 

^ -2 -1 fl~ T~ 2" ^ -2~ -1 Y 2 

10 10 10 10 10 10 10 10 10 10 10 10 
t/X 

FIG. 4: Non-Gaussian parameter 0:2(0 for the layer-to layer (#) and in-layer (O) diffusion for 
smectic LCs of hard spherocylinders with varying length-to-diameter ratios L* and reduced pres- 
sures P*. Length-to-diameter ratio, pressure, and packing fraction are given in each frame as (L*, 
P*, r?). 

Most of the information obtained by the analysis of the non-Gaussian parameter can also 
be deducted by the mean square displacements (MSDs), (A;^^^^)) and {Ax^ (t) + Ay^ (t)) , 
shown in Fig. 5. The xy-MSD (open circles) is characterized by a relatively smooth crossover 
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FIG. 5: Mean square displacement (MSD) in units of along the nematic director n (solid circles) 
and along the plane perpendicular to n (open circles) for smectic LCs of hard sphero cylinders with 
varying length-to-diameter ratios L* and reduced pressures P* . Length-to-diameter ratio, pressure, 
and packing fraction are given in each frame as (L*, P*, r/). The solid lines are a guide for the eye. 

from the short- to long-time diffusion, as observed in slightly dense liquids. By contrast, for 
the 2:-MSD (solid circles) one clearly detects a more sophisticated behavior as three separate 
time regimes can be identified. The short-time dynamics, with the rods still rattling in 
their cages, are diffusive, that is (Az^(t)) oc t. In this regime, {^z^{t)) > (Ax'^it) + Ay'^{t)) 
because of the anisotropy of the rods [56j . After an induction time, we observe the formation 
of a plateau which extends up to t^ax and quantifies the time to escape from the trapping 
cages. In this time window, the dynamics becomes subdiffusive. Finally, at different times, 
the xy-MSD and z-MSD become linear with time and the long-time diffusive regime is 
reached. 

From the MSDs in the diffusive regime, we computed the long-time diffusion coefficients in 
the z- and x?/-directions by applying the well-known Einstein relation [4j. The values of the 
long-time diffusion coefficients, D^^ and D^, are presented in Table I in units of D'^/t. The 
dynamics of each system is characterized by a diffusion coefficient in the 2;?/-direction that 
is larger than the one in the z-direction. This result is in agreement with the dynamics in 
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0.56 



FIG. 6: Diffusion ratio 70 (top) and dimensionless diffusion coefficients (bottom) as a function of 
the packing fraction for smectic LCs of liard splierocylinders witli length-to-diameter ratio L* = 5.0. 
^*xy,z — ^xy,z'^/^'^ denotes the reduced diffusion coefficients, as given in Table I. The solid lines 
are power law fits. 

thermotropic smectogenic LCs [US], but in contrast with recent experiments on the diffusion 
of fd viruses [2Z], most probably because of their huge aspect ratio (L* > 100). In Fig. 6, 
we exemplarily show the dependence on the packing fraction of the diffusion coefficients and 
diffusion ratio, defined as 



7o 



xyl 



(6) 



for L* = 5.0. D^y and as well as 70 are well fitted by power law functions of the type 
1]'", with u ~ 4.6, 26.0, and -21.4, respectively. Increasing the pressure has a significantly 
larger effect on the inter-layer dynamics than on the in-layer one, as decreases much 
faster than -D^y This is to be expected as the energy barriers (see Fig. 2) hamper the 
diffusion more in dense systems. For the lowest packing fractions, i.e., where the smectic 
phase almost coexists with the nematic phase, the two diffusion coefficients approach each 
other. Similar considerations are also valid for the systems with aspect ratio L* = 3.4 and 
3.8 (not shown here). 

In addition, we analyzed the structural relaxation by calculating the self-part of the 
intermediate scattering functions, defined in Eq. (4). In Fig. 7, we show Fg^xyit) and 
Fs,zif) at the wave vectors q = (gx,Q'y,0) and (0,0,^2), respectively, corresponding to the 
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main peaks of the static structure factor. We found that {q^ + Qy) ~ 6 for all systems, 
whereas Dq^ — 1.4, 1.3, and 1.0 for the systems with L* = 3.4,3.8, and 5.0, respectively. 
Regardless the aspect ratio, we can affirm that the in-layer structural relaxation is several 
orders of magnitude faster than the inter-layer relaxation. If we define the relaxation time 
as the time at which Fs{t) decays to e~^, then tr^xy/T is of the order of 10~^ and tr^z/'T > 10^. 
We also find that Fs^xy{t) decays very fast to zero with slightly stretched exponential decay, 
as expected for dense liquid-like dynamics [66] . 

By contrast, the inter-layer relaxation develops in two steps separated by a plateau, the 
beginning of which corresponds to the development of the cage regime. During the initial 
decay, which is relatively fast {t/r < 1), the rods are free to rattle inside the temporary cage 
formed by the nearest neighbors [38j, without perceiving their trapping effect, and Fs^z{t) is 
characterized by an exponential decay. After this short time lapse, we detect a plateau, whose 
height and temporal extension depend on the packing of the system, as was also found in 
colloidal glasses [66]. In our previous analysis, we observed that the Gaussian approximation 
of the self-intermediate scattering function, that is Ff^{t) = exp[—q'^{Az'^(t))], does not 
show any significant plateau |45j. This result indicates that the existence of a plateau 
must be linked to the non- vanishing NGP 02,2, and hence to the heterogeneous inter- layer 
dynamics. After the plateau, a second decay, which corresponds to the escape from the 
temporal cages, leads the systems towards the structural relaxation on a time scale which is 
long at the highest packing fractions. We were only able to estimate the long-time relaxation 
decay in the z-direction for systems 5*3 and S^, due to their relatively low packing fractions. 
In particular, we fit the long-time decay of Fg^zit) with a stretched exponential function of 
the form exp [-{t/U)^] , with U/t = 2500 and /3 = 0.8 for ^3, and U/t = 650 and /3 = 0.6 
for 5*5. The relaxation time of the remaining systems are beyond our simulation time. 

The results shown so far give clear evidence of the existence of fast-moving particles 
determining a distribution of decay rates which affects the long-time structural relaxation. 
We now turn our attention to the occurrence of collectively moving particles which might 
play a crucial role, or might even be responsible, for the heterogeneous inter-layer dynamics 
in smectic LCs. To this end, we ffist label the fast-moving particles by applying the static 
cluster criterion defined below Eq. (5). Fig. 8 shows the probability size distribution P{n) 
of the number of interlayer rods n in a string-like cluster using rattling distance Srat = 2o", 
where a is the standard deviation as specified in Eq. (5) and presented in Table I. In 
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FIG. 7: Self-intermediate scattering function for smectic LCs of hard spherocylinders with length- 
to-diameter ratio L* = 3.4, and reduced pressures P* = 2.85 and 3.00 (top), L* = 3.8, and 
P* = 2.35 and 2.50 (middle), and L* = 5.0, and P* = 1.60 and 2.00 (bottom). The solid and 
open symbols refer to the lowest and highest pressure, respectively. Squares and circles refer, 
respectively, to the in-layer and inter-layer relaxation. The solid lines are fits. 

Fig. 9, we give an illustrative example of static strings observed in system 5*3. Regardless 
the particle anisotropy, two interesting conclusions can be drawn: (z) the observed strings 
consist mostly of 2 or 3 rods, while clusters of more than 5 rods are extremely rare, but 
do exist; and iii) strings containing more than 2 rods are more often formed in the denser 
state, as usually observed in supercooled liquids and glassy systems, where the average 
cluster size increases when the caging effect becomes stronger [IH]. The semilog plot of 
Fig. 8 shows that the size probability distribution is roughly exponential. In particular, 
we observed that Pin) oc exp (— ari), from which the average cluster size can be estimated: 
(n) = (1 — exp(— a))^^. The largest clusters are thus found for the shortest rods and the 
highest pressure. 

The energy barriers of Fig. 2 and the periodically peaked shape of the VHFs of Fig. 3 
unfold the rattling-and-jumping layer-to-layer diffusion of the rods, but only provide a global 
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FIG. 8: Size probability distribution P(n) of the number of interlayer rods n in a stringlike clusters 
(with 5rat = 2o", see text), in a smectic LCs of hard spherocylinders with length-to-diameter ratio 
L* = 3.4, and reduced pressures P* = 2.85 and 3.00 (top), L* = 3.8, and P* = 2.35 and 2.50 
(middle), and L* = 5.0, and P* = 1.60 and 2.00 (bottom). The solid and open symbols refer to the 
lowest and highest pressure, respectively. The solid lines denote the fit P{n) oc exp(— an), with a 
given in the figure. 

picture of the dynamics in smectic LCs. In order to gain a deeper understanding of the actual 
dynamics on the particle scale behind the layer-to-layer hopping-type diffusion, we followed 
the trajectories of single particles. Interestingly, we observed that some rods diffuse very 
fast, others move to the inter-layer spacing, where they might reside for a long time, and 
then return to their original layer, and others move from one layer to another several times 
or perform consecutive jumps as shown in Fig. 10b for system S^. Furthermore, transverse 
inter- layer rods, although extremely rare due to the high energy barriers [GT]. [62] , have also 
been detected. These transverse inter-layer particles might diffuse either to a new layer or 
go back to the old one by keeping or changing their original orientation. This variegated 
behavior suggests a rather broad distribution of layer-to-layer jump times n(tj), where tj 
is the time between the first and last contact with the new and old layer, respectively. Such 
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FIG. 9: (color online). Snapshot of 3000 rods with length-to-diameter ratio L* = 3.8 and packing 
fraction ij = 0.536. In-layer rods and single interlayer rods (black), both with diameters reduced 
to I?/4 for clarity, are predominantly shown. The thicker rods denote transverse ones (red or dark 
shaded) as well as stringlike clusters of 2 (yellow or light shaded) up to 11 (brown or dark shaded) 
rods. 

a contact is established as soon as the particle is at a rattling distance 6rat = 2cr from the 
middle of a smectic layer. The probability distributions of jump times for systems Si — Sq 
have been computed by averaging over at least N jumps, with N the number of rods in each 
system. In Fig. 10a, we give the probability distribution of jump times n(tj) for system 
5*3 . As expected, the distribution is not particularly narrow, but extends over two time 
decades 0.01 < tj/r < 1, with the most probable jump time ^=0.17r and the median time 
(i.e. the time at which 50% of the jumps have been performed) t*j = 0.28r. These times 
increase at rj = 0.551, as given in Table I. The remaining systems with shorter or longer 
rods show the same tendency and similar times (see Table I). We also distinguished single 
and consecutive (or multiple) jumps according to the dwelling time being longer or shorter, 
respectively, than t*j before another jump is started. Single jumps are significantly more 
frequent than multiple jumps, especially at high densities where the latter are less than 1% 
of the total number of jumps. Most of the multiple jumps consist of double or triple jumps, 
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FIG. 10: (color online), (a) Probability distribution of jump times n(tj) based on 5rat = 2(T for 
lenght-to-diameter ratio L* = 3.8 and packing fraction rj = 0.536. (b) Trajectories of jumping rods 
in the same system projected onto the xz plane, with the dashed lines locating the centre of the 
smectic layers. The arrow indicates the trajectory of a transverse interlayer particle. 

while quadruple ones are basically irrelevant. The largest number of multiple jumps was 
observed in the system with the lowest packing fraction, i.e., : L* = 5.0,1] = 0.508, where 
~ 5% of the total jumps is multiple, of which ^ 83% is a double and the remaining fraction 
is a triple jump. 

The formation of static stringlike clusters does not necessarily imply the occurrence of 
dynamic cooperativity, as a rod that belongs to a cluster might still diffuse individually from 
layer-to-layer or might fail to jump. In order to ascertain the occurrence of collective motion 
of strings, we introduced a dynamic cluster criterion. More specifically, we assume that two 
jumping rods i and j are actually moving cooperatively, if (1) their arrival times t*^*-* and t^^^ 
in their new layers satisfy the condition \ f — P\ < Ato and (2) their z and xy distances at t^*^ 
and t(^) are smaller than h and D, respectively. To select a consistent value for Ato, we use 
the distribution of jump times n(tj). If we assume that Atg is the maximal jump time t'J'"'^ 
(i.e. long enough for all jumps to be performed according to n(tj)), we find that for L* = 3.4 
the ratio between the number of collective jumps and the total number of jumps is as 
given in Table I. These jumps involve mainly two or three rods, whereas collective jumps of 
4 or more rods are extremely rare. We further note that the vacated space of a jumping rod 
can be either occupied by another rod jumping in the same direction or by a rod jumping in 
the opposite direction, with roughly the same probability. A less restrictive spatial criterion 
would not affect significantly these values. By contrast, the number of collective jumps is 
rather sensitive to the temporal criterion. If we reduce Ato to the median jump time t*j, then 
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the fraction of collective jumps, /*, decreases substantially as shown in Table I. Regardless 
the details of the cluster criterion, we thus find a fraction of 10^^ — 10^^ of the jumps to be 
collective, the more so for longer rods at lower pressures. 

The motion is indeed strongly cooperative at low packing fractions, despite the larger 
static stringlike clusters detected in the denser systems (see Fig. 8). This result is most prob- 
ably due to the permanent smectic barriers which increase upon approaching the smectic-to- 
solid phase transition and hence hamper the attempted jumps of the rods in the strings. By 
contrast, in glass-forming systems, where no permanent barriers are observed, the cluster 
size increases upon approaching the glass transition EZ] • 

IV. CONCLUSIONS 

In summary, we have studied the diffusion and structural relaxation in equilibrium smec- 
tic LC phases of hard rods with different anisotropics by computer simulations. Remarkably, 
these systems exhibit non-Gaussian layer-to-layer diffusion and dynamical heterogeneities 
which are similar to those observed in out- of- equilibrium supercooled liquids. The simulta- 
neous presence of temporary cages due to the trapping action of neighboring rods and the 
permanent barriers due to the static smectic background, provokes a rattling-and-jumping 
diffusion which influences the long-time structural relaxation decay. In analogy with glassy 
systems, one can clearly distinguish three separate time regimes for z-motion. The short- 
time diffusive regime, with the rods rattling around their original location without feeling the 
surrounding neighbors, is characterized by a Gaussian distribution of the z-displacements 
and an exponential temporal relaxation. The subdiffusive regime at intermediate times dis- 
plays a non-Gaussianity and a plateau in both mean square displacement and intermediate 
scattering function. At this stage, the interlayer dynamics is heterogeneous with fast-moving 
particles diffusing individually or cooperatively in a stringlike fashion. Finally, at long times, 
the systems enter a second diffusive regime with Gaussian distributions of the displacements 
and non-exponential decay of the intermediate scattering function. By contrast, at all time- 
scales, the in-layer diffusion is typical for two-dimensional fluids with a negligible NGP and 
a structural relaxation which is at least 4 time decades faster than the inter-layer one. 

The analysis of the self-VHFs points out the tendency for the rods to diffuse from layer 
to layer through quasi-discretized jumps. This hopping-type motion is significantly ham- 
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pered in very dense systems, where the barriers for the layer-to-layer diffusion intensify the 
confinement of the rods to the middle of the smectic layers. The temporal extension of the 
jumps is not uniform, but characterized by a rather broad time distribution which covers 
approximately two time decades. Depending on the dwelling time between two successive 
jumps of the same particle, single and multiple jumps have been detected, with the former 
significantly more frequent than the latter, especially at high densities. Although the inter- 
layer rods are usually oriented along the nematic director, some of them assume a transverse 
orientation. The long tails of the VHFs indicate the presence of particles that are able to 
diffuse much longer distances than the average, especially at low packing fractions. Inter- 
estingly, the dynamic behavior of such fast-moving particles supports the intriguing analogy 
with glassy systems even further. In particular, fast-moving particles assemble in string- 
like clusters whose average length increases upon approaching the smectic-to-crystal phase 
transition. Likewise, fragile and strong glass-formers show a similar tendency when cooled 
down towards the glass transition temperature [ISl [32] • We also gave clear evidence that 
the strings detected in static configurations can promote collective diffusion of jumping par- 
ticles. This is especially tangible at low packing fractions where the hampering action of the 
permanent energy barriers is less effective. Finally, we also investigated the effect of particle 
anisotropy on the non-Gaussian layer-to-layer diffusion and cooperative motion in smectic 
LCs. We find that at fixed packing fraction, the barrier height increases with increasing par- 
ticle anisotropy, and hence the dynamics is more heterogeneous and non-Gaussian for longer 
rods, yielding a lower diffusion coefficient along the nematic director and smaller clusters 
of inter-layer particles that move less cooperatively. At fixed barrier height, the dynamics 
becomes more non-Gaussian and heterogeneous for longer rods; smaller clusters move more 
collectively, giving rise to a higher diffusion coefficient along the nematic director. 

Our results, which might be relevant for the study of the dynamics in confined fluids 
[55] as well as for the diffusion of lipids and proteins in cellular membranes p5| [70] . are 
already stimulating the analysis of dynamical processes in columnar liquid crystals, where 
the presence of inter-columnar energy barriers provokes a non-Gaussian hopping-type mo- 
tion and a relaxation behavior that is also remarkably similar to that of out-of-equilibrium 
supercooled liquids |7lJ. We hope that our findings stimulate further theoretical and exper- 
imental studies of particle-scale dynamics in heterogeneous (confined, smectic, columnar) 
liquids. 
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